Homogenization and asymptotics for small transaction 

costs 

H. Mete S ONER* Nizar Touzi''' 

March 16,2012 

Abstract 

We consider the classical Merton problem of lifetime consumption-portfolio optimiza- 
tion problem with small proportional transaction costs. The first order term in the asymptotic 
expansion is explicitly calculated through a singular ergodic control problem which can be 
solved in closed form in the one-dimensional case. Unlike the existing literature, we consider 
a general utility function and general dynamics for the underlying assets. Our arguments are 
based on ideas from the homogenization theory and use the convergence tools from the the- 
ory of viscosity solutions. The multidimensional case is studied in our accompanying paper 
1311 using the same approach. 
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1 Introduction 



The problem of investment and consumption in a market with transaction costs was first studied 
by Magill & Constantinides 1,26] and later by Constantinides jlO|. Since then, starting with 
the classical paper of Davis & Norman ifTTI an impressive understanding of this problem has 
been achieved. In these papers and in lfT2l[36l the dynamic programming approach in one space 
dimension has been developed. The problem of proportional transaction costs is a special case 
of a singular stochastic control problem in which the state process can have controlled discon- 
tinuities. The related partial differential equation for this class of optimal control problems 
is a quasi-variational inequality which contains a gradient constraint. Technically, the multi- 
dimensional setting presents intriguing free boundary problems and the only regularity result to 
date are 1341 and fJSl . For the financial problem, we refer to the recent book by Kabanov & 
Safarian |[22l . It provides an excellent exposition to the later developments and the solutions in 
multi-dimensions . 

It is well known that in practice the proportional transaction costs are small and in the 
limiting case of zero costs, one recovers the classical problem of Merton [28]. Then, a natural 
approach to simplify the problem is to obtain an asymptotic expansion in terms of the small 
transaction costs. This was initiated in the pioneering paper of Constantinides ifTOl . The first 
proof in this direction was obtained in the appendix of li36J . Later several rigorous results 
in [18] |20l |32l and formal asymptotic results jl] [19] [38ll have been obtained. The rigorous 
results have been restricted to one space dimensions with the exception of the recent manuscript 
by Bichuch and Shreve [6]. 

In this and its accompanying paper ifSTl . we consider this classical problem of small propor- 
tional transaction costs and develop a unified approach to the problem of asymptotic analysis. 
We also relate the first order asymptotic expansion in e to an ergodic singular control problem. 

Although our formal derivation in Section[3]and the analysis of 1311 are multi-dimensional, 
to simplify the presentation, in this introduction we restrict ourselves to a single risky asset 
with a price process {St,t > 0}. We assume St is given by a time homogeneous stochastic 
differential equation together with So = s and volatility function (j(-). For an initial capital 
z, the value function of the Merton infinite horizon optimal consumption-portfolio problem 
(with zero-transaction costs) is denoted by v{s,z)- On the other hand, the value function for 
the problem with transaction costs is a function of s and the pair {x,y) representing the wealth 
in the saving accounts and in the stock. Then, the total wealth is simply given by z — x + y. 
For a small proportional transaction cost > 0, we let v^{s,x,y) be the maximum expected 
discounted utility from consumption. It is clear that v^{s,x,y) converges to v{s,x+y) as e tends 
to zero. Our main analytical objective is to obtain an expansion for v*^ in the small parameter e. 

To achieve such an expansion, we assume that v is smooth and let 



be the corresponding risk tolerance. The solution of the Merton problem also provides us an op- 
timal feedback portfolio strategy 9{s,z) and an optimal feedback consumption function c{s,z)- 
Then, the first term in the asymptotic expansion is given through an ergodic singular control 
problem defined for every fixed point {s,z) by 



a{s,z) := infy(,s',z,M), 



M 



where M is a control process of bounded variation with variation norm \\M\ 



J(s,z,M) := limsup [ M^llit + \\M\\r 



T^oo 1 [Jo 
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and the controlled process ^ satisfies the dynamics driven by a Brownian motion B, and param- 
eterized by the fixed data {s,z): 



d^t [B{s,z){l - d-{s,z))]dB, + dM,. 



The above problem is defined more generally in Remark [331 and solved explicitly in the sub- 
section |4T|below in terms of the zero-transaction cost value function v. 

Let {Z^" , f > 0} be the optimal wealth process using the feedback strategies Q , c, and starting 
from the initial conditions Sq = s and Zq ~ = z- Our main result is on the convergence of the 
function 



Main Theorem. Let a be as above and set a := r\v,a. Then, as e tends to zero, 



u^{x,y) — ^ u{s,z) := E 



e-P'a{S,,Z''')dt 



locally uniformly. 



(1.2) 



Naturally, the above result requires assumptions and we refer the reader to Theorem l6.1l for 
a precise statement. Moreover, the definition and the convergence of is equivalent to the 
expansion 

v'{s,x,y)=viz)~e\,{z) + o{e-), (1.3) 

where as before z=x+y and o(£*^) is any function such that o[e'^)/e'^ converges to zero locally 
uniformly. 

A formal multi-dimensional derivation of this result is provided in Section[3] Our approach 
is similar to all formal studies starting from the initial paper by Whalley & Willmont lf38l . 
These formal calculations also provide the connection with another important class of asymp- 
totic problems, namely homogenization. Indeed, the dynamic programming equation of the 
ergodic problem described above is the corrector equation in the homogenization terminology. 
This identification allows us to construct a rigorous proof similar to the ones in homogenization. 
These assertions are formulated into a formal theorem at the end of Section [3] The analysis of 
Section|3]is very general and can easily extend to other similar problems. Moreover, the above 
ergodic problem is a singular one and we show in lUTj that its continuation region also describes 
the asymptotic shape of the no-trade region in the transaction cost problem. 

The connection between homogenization and asymptotic problems in finance has already 
played an important role in several other problems. Fouque, Papanicolaou & Sircar ll24l use 
this approach for stochastic volatility models. We refer to the recent book jZSl for information 
on this problem and also extensions to multi dimensions. In the stochastic volatility context the 
homogenizing (or the so-called fast variable) is the volatility and is given exogenously. Indeed, 
for homogenization problems, the fast variable is almost always given. In the transaction cost 
problem, however, this is not the case and the main difficulty is to identify the "fast" variable. 
A similar difficulty is also apparent in a problem with an illiquid model which becomes asymp- 
totically liquid. The expansions for that problem was obtained in ||30ll . We use their techniques 
in an essential way. 

The later sections of the paper are concerned with the rigorous proof. The main technique is 
the viscosity approach of Evans to homogenization lfT3l[T4l . This powerful method combined 
with the relaxed limits of Barles & Perthame 12J provides the necessary tools. As well known, 
this approach has the advantage of using only a simple U° bound which is described in Section 
12] In addition to 121 [131 11411 . the rigorous proof utilizes several other techniques from the theory 
of viscosity solutions developed in the papers ||2l [15| [l2l|23] |33] [37l for asymptotic analysis. 
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For the rigorous proof, we concentrate on the simpler one dimensional setting. This simpler 
setting allows us to highlight the technique with the least possible technicalities. The more 
general multi-dimensional problem is considered in ifsTl . 

The paper is organized as follows. The problem is introduced in the next section and the 
approach is formally introduced in Section[3] In one dimension, the corrector equation is solved 
in the next section. We state the general assumptions in Section |5] and prove the convergence 
result in Section |6l In Section [T] we discuss the assumptions. Finally a short summary for the 
power utility is given in the final Section. 



2 The general setting 

The structure we assume is the one developed and studied in the recent book by Kabanov & 
Safarian 1221 . We briefly recall it here. 

We assume a financial market consisting of a non-risky asset 5''^ and d risky assets with price 
process {St ~ {S} . . ,Sf),t > 0} given by the stochastic differential equation, 

-g- = r(St)dt, 
dS' 

-f = ^'iS,)dt + Y,a'-JiS,)dW/, l<i<d, 
\ j=i 

where r : R'' i-^ IR+ is the instantaneous interest rate function, and 

fi-.W^^W^, and a ■.Mf' ^ 

are the coefficients of instantaneous mean return and volatility. We use the notation ^j(R) 
to denote d x d matrices with real entries. Under the standard Lipschitz and linear growth 
conditions (which we assume without stating) the above stochastic differential equation has a 
unique strong solution. 

The portfolio of an investor is represented by the dollar value X invested in the non-risky 
asset and the vector process Y = {Y\... of the value of the positions in each risky asset. 
The portfolio position is allowed to change in continuous-time by transfers from any asset to 
any other one. However, such transfers are subject to proportional transaction costs. 

We continue by describing the portfolio rebalancing in the present setting. For all i,j = 
0,...,d, let L','^ be the total amount of transfers (in dollars) from the i-th to the j-th asset cumu- 
lated up to time t. Naturally, the processes {L','\t > 0} are defined as cad-lag, nondecreasing, 
adapted processes with Lq- ~ and U ' = 0. The proportional transaction cost induced by a 
transfer from the i-th to the j-th stock is given by e^X''^ where e > is a small parameter, and 

X'j>0, A''=0, ij = 0,...,d. 

The scaling £^ is chosen to state the expansion results simpler We refer the reader to the recent 
book of Kabanov & Safarian 1221 for a thorough discussion of the model. 

The solvency region Kg is defined as the set of all portfolio positions which can be trans- 
ferred into portfolio positions with nonnegative entries through an appropriate portfolio rebal- 
ancing. We use the notation {(''■')ij=o,...d G -^d+i ) to denote this appropriate instantaneous 
transfers of size £'■■'. Then, Kg is given by. 

Kg := I {x,y) eRxR': {x,y) + E ( {i^-' - (1 + e'X''^)e''^)) £ R\+' 

for some (f •■')o<ij<d e ^d+i(M+)|, 
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where {eo,. . . ,6^) denotes the canonical basis of M''+'. It is clear that for any it would not 
be optimal to have both £''j > and f-' ' > 0. For later use, we also define 

Afj -.^ei-ej + e^X'-' ei, ij^0,...,d. 

In addition to the trading activity, the investor consumes at a rate determined by a nonnega- 
tive progressively measurable process {c,,r > 0}. Here c, represents the rate of consumption in 
terms of the non-risky asset S^. Such a pair v := (c,L) is called a consumption-investment strat- 
egy. For any initial position {Xq-,Yq-) = {x,y) G R x M.'', the portfoUo position of the investor 
are given by the following state equation, 

d 

dXt = {r{S,)Xt - c,)dt + 52 [dLi-^ - (1 + £^X''^-^)dl}l^ ) , 

./=i 

j=Q ^ ^ 

The above solution depends on the initial condition {x,y), the control v and also on the initial 
condition of the stock process s. Let (X,}')^-* '^'^ be the solution of the above equation. Then, a 
consumption-investment strategy v is said to be admissible for the initial position {s,x,y), if 

{X,Y)J-''-^' eKs, Vf>0, P-a.s. 

The set of admissible strategies is denoted by &^{s,x,y). For given initial positions Sq — sE R'^, 
Xq- = X e M, Yq- = y e M'', the investment-consumption problem is the following maximization 
problem, 



v^{s,x,y) := sup E 

{c.L)e&Hs.x.,y) 



U{c,)dt 



where C/ : (0,°°) ^ M is a utility function. We assume that U is C^, increasing, strictly concave, 
and we denote its convex conjugate by, 

t7(c) sup {f/(c)-cc} , cGM. 

OO 

Then f7 is a convex function. It is well known that the value function is a viscosity solution 
of the corresponding dynamic programming equation. In one dimension, this is first proved in 
ll36l . In the above generality, we refer to ll22l . To state the equation, we first need to introduce 
some more notations. We define the second order linear partial differential operator, 

^ := M-(D.v + D,0+'-Dx + ^Tr[aaT(D,,, + D,, + 2D„)], (2.1) 

where ^ denotes the transpose and for /, y = I,. . . ,d, 

D..= .-, D,.= .^, Dv.= y^, 

^2 . . ^2 ■ ■ d- 

ds'dsJ' » dy'dyJ' ^ ds'dyJ' 

D.V = (Ds)l<,<rf, Dy = (Dy)l<Krf, Dyy := (Dyy )l<ij<rf, Djs := (Df/ ) l<ij<rf, D.vy := {T>sy )\<iJ<d- 

Moreover, for a smooth scalar function {s,x,y) £ R'{_ x R x R'^ i— >■ (p{x,y), we denote: 

(p,:=^ eR, (p.-=^ eR''. 

ox ay 
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Theorem 2.1 Assutne that the value function is locally bounded. Then is a viscosity 
solution of the dynamic programming equation in x Ks, 

, min { j3 - ^v' - U{vl) , Af.,- • (v^, v,^) } = 0. (2.2) 

Moreover, v^ is concave in {x,y) and converges to the Merton value function v := v^, as e > 
tends to zero. 

Under further conditions the uniqueness in the above statement is proved in ll22ll . However, 
this is not needed in our subsequent analysis. 



2.1 Merton Problem 

The limiting case of £ = corresponds to the classical Merton portfolio-investment problem in 
a frictionless financial market. In this limit, since the transfers from one asset to the other are 
costless, the value of the portfolio can be measured in terms of the nonrisky asset 5". We then 
denote hy Z: = X + + . . .+Y^' the total wealth obtained by the aggregation of the positions on 
all assets. In the present setting, we denote by d' :=Y' and d := (0' , . . . , 0'') the vector process 
representing the positions on the risky assets. The wealth equation for the Merton problem is 
then given by 

dZt = (r(5,)Z, - cr)dt + L ^/ V - • (2.3) 

An admissible consumption-investment strategy is now defined as a pair (c, 0) of progressively 
measurable processes with values in M+ and , respectively, and such that the corresponding 
wealth process is well-defined and almost surely non-negative for all times. The set of all 
admissible consumption-investment strategies is denoted by 0(i,z). 
The Merton optimal consumption-investment problem is defined by 



v(s,z) '.— sup E 

(c,e)G0(s„-:) 



U{c,)dt 



Throughout this paper, we assume that the Merton value function v is strictly concave in z and 
is a classical solution of the dynamic programming equation. 



J3y-rzy,- Jf°v-f/(v-)- sup • ((jU - rlrf)v, + (7(7^D,,v) + i|(7^0pv-,| = 0, 

em'' ^ 

where 1</ := (1, . . . , 1) e W' , D,, |:D,, and 



A'-D. + ^Tr[(J(7TDj. (2.4) 

The optimal controls are smooth functions c{s,z) and d{s,z) obtained by the maximization of 
the Hamiltonian. Hence, 

= j3v-^°v-f7(v,) (2.5) 

-rzvz -e-in- rld)v, - aa'^e ■ D,-v - ^ | a'^e\^v„, 

the optimal consumption rate is given by, 

c{s,z) := -U'{vzis,z)) ^{U'y\v,{s,z)) for seE.i,z>0, (2.6) 
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and the optimal investment strategy 9 is obtained by solving the finite-dimensional maximiza- 
tion problem, 

max lha'^9\^v-, + • ((jU - r)v- + c7(7'^D,,v) |. 
Since v is strictly concave, the Merton optimal investment strategy 9{s,z) satisfies 

-v„{s,z) a{s)cy'^{s)9{s,z) = {fi - rld){s)v,{s,z) + a{s)a'^D,,v{s,z). (2.7) 

3 Formal Asymptotics 

In this section, we provide the formal derivation of the expansion in any space dimensions. 
In the subsequent sections we prove this expansion rigorously for the one dimensional case. 
Convergence proof in higher dimensions is carried out in a forthcoming paper ||3T|. In the 
sequel we use the standard notation 0{e^) to denote any function which is less than a locally 
bounded function times and o(e'-) is a function such that o(e'^)/e* converges to zero locally 
uniformly. 

Based on previous results IMlfflllllEQlElllMl, we postulate the following expansion, 

v''{s,x,y) = v{s,z) - e^u{s,z) - eV(i,z,^) + o{e-), (3.1) 
where {z,£,) — {z,£,£) is a linear transformation of {x,y) G Kg given by 

z=x + / + ...+y", <^':=^;(x,y) = ^:— /=!,.. 

= (0 ' , . . . , 0'') is the Merton optimal investment strategy of i2.7\ . In the postulated expansion 
(13. Il l, we have also introduced two functions 

M : R'l X R+ M, and w : M'^ x R+ x R'' ^ R. 

Notice that the above expansion is assumed to hold up to e^, i.e. the o[e^) term. Therefore, 
the reason for having an higher term like e^w{z, t, ) explicitly in the expansion may not be clear 
However, this term contains the fast variable t, and its second derivative is of order e^. This 
follows the intuition introduced in the pioneering work of Papanicolaou and Varadhan 1291 in 
the theory of homogenization. 

The main goal of this section is to formally derive equations for these two functions. A 
rigorous proof will be also provided in the subsequent sections and the precise statement for 
this expansion is stated in Section |6l 

Since {x^y) G /Tg i— >• (z,^) G R+ x M'' is a one-to-one change of variables, in the sequel for 
any function / of {s,x,y) we use the convention, 

f{s,z,^):=f{s,z^e^-9{s,z),£^+e{s,z)). (3.2) 

The new variable ^ is the "fast" variable and in the limit it homogenizes to yield the convergence 
of v^(s,z,^) to the Merton function v{s,z) which depends only on the (5,z)-variables. This is 
the main formal connection of this problem to the theory of homogenization. This variable 
was also used centrally by Goodman & Ostrov |19j. Indeed, their asymptotic results use the 
properties of the stochastic equation satisfied by e^'^(X,,y,). 

We continue by formally deriving the equations satisfied by u and w. First we directly 
differentiate the expansion (13. It and compute the terms appearing in ( 12. 2t in term of u and w. 
The directional derivatives are given by. 
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We directly calculate that, 

{wx{s,z,^),Wy{s,z,^)) = (w,~^0,-w^^ld+i + j{0,w^) . (3.3) 
To simplify the notation, we introduce 

D^wis,zA) {0,D^wis,z,^)) e R''+^ (3.4) 

Then, 

Kj ■ i^'^^y) = (^'''^'^ + ('J - ■ + 0{£')- (3.5) 

The elliptic equation in ( 12.2b requires a longer calculation and we will later use the Merton 
identities ( |231 ). ( |Z6] l and ^J}. Firstly, by ( |23T l. 

= (0-y)- [(Ai-rl^)v, + CTCTX:V'] +^(|(tT0|2-|(tT3;|2)v,, 

-E^(pu-^u) + ^Ti[aa'^T)yyw] + 0{e^). 



We use a Taylor expansion on the terms involving U. and ( |2.6l )-( |Z7] i in the first line. The result 
is, 



£4 



-e- ( jSw - + cii, ) + y Tr[(JC7^Dvvw] + 0{e^) 



1 / \ e 

/I \ 

E^i^--\a'^^\\,-l5u + ^u-cu,j + —Ti[aa'^T)yyw] + 0{e^). (3.6) 



Finally, from ( 13.31 ). we see that 

1 



dyW = w.lj + -{ld- lf/0r ) Wt, ■ 



Therefore, 



We substitute this in ( I3.6l l and use the fact that y = + (9(e). This yields, 

I"" ^£^[-\\(y^^?-v-z + ]^Tx[aa^w^^] -s^u^+0{e^), (3.7) 
where a{s,z) is given by 

ais,z) = {Id - 0,(*,z)i;^)diag[0(^,z)](j(.), (3.8) 
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diag[0] denotes the diagonal matrix with i-th diagonal entry 9', and 

s!/u = Pu- -{rz+e-{ii- r\,i) - c)u- - ^ P "zz - c^c^^^ ' Dk"- (3.9) 

Observe that the above operator is the infinitesimal generator of the pair process (5,2) where Z 
is the optimal wealth process in the Merton zero-transaction cost problem corresponding to the 
optimal feedback controls (c, 0). In particular, the dynamic programming equation ( 12.51 ) for the 
Merton problem may be expressed as, 

£!^v{s,z)=U{c{s,z)). (3.10) 

We have now obtained expressions for all the terms in the dynamic programming equation ( 12.21 ). 
We substitute ( I3.5l l and ( I3.7l i into ( 12.2b . The result is the following equation for w and u. 



max max \^l-\o^ {s)B,fv,,{s,z) ~\Tr[aa^ {s,z)w^^{s,z,B,)\ + a{s,z) , 

-A''V,(^,z) + (e,-e,).D^w(z,^)} = 0. 
where = (0,D^w) is as in ( 13.41 ) and a is simply defined by, 

fl(s,z) :=£/M(i,z), ^eR'l, z>0. 

In the first equation above, the pair (5,2) is simply a parameter and the independent variable is 
^. Also the value of the function w(i,z,0) is irrelevant in ( 13.11 ) as it only contributes to the 
term. Therefore, to obtain a unique w, we set its value at the origin to zero. We continue by 
presenting these equations in a form that is compatible with the power case. So we first divide 
the above equation by and then introduce the new variable p —^lr\(s^z) where 77 is the risk 
tolerance coefficient defined by ( 11.11 ). We also set 

w{s,z,p)~ ' N / N ' ^ ^'^ \ t "O^'^ :=rW- 

77(i,z)v,(i,z) ri{s,z)vz{s,z) ri(s,z) 

Then, the corrector equations in this context is the following pair of equations. 



Definition 3.1 (Corrector Equations) For a given point {s,z) G M.'[ x R^, the first corrector 
equation is for the unknown pair (a(i,z),w'(i,z, •)) G R x C^iW'), 

max max<^-J - -Tr[aci^ {s,z)wpp{s,z,p)\ + a{s,z) , (3.11) 

Q<i.j<d y L L 

-A'-' + (e,-e^)-DpM>(i,z,p)} =0, VpeM'', 

together with the normalization h>(s,z,0) = 0. 

The second corrector equation uses the constant term fl(s,z) from the first corrector equation 
and it is a simple linear equation for the function 

^M(i,z) =a(i,z) = v,(i,z)77(i,z)fl(i,z), VieR+, zeR+. (3.12) 

We say that the pair is the solution of the corrector equations for a given utiUty function 
or equivalently for a given Merton value function. FDJ 
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We summarize our formal calculations in the following. 



Formal Expansion Theorem. The value function has the expansion ( 13.1b where is the 

unique solution of the corrector equations. 

Remark 3.1 The function u introduced in ( |1.2t is a solution of the second corrector equation 
( 13.121 ). provided that it is finite. Then, assuming that uniqueness holds for the linear PDE ( |3.12t 
in a convenient class, it follows that u is given by the stochastic representation ( 11.2b . 

Remark 3.2 Usually a second order equation like ( 13.121 ) in (0,oo) needs to be completed by a 
boundary condition at the origin. However, as we have already remarked, the operator jz/ is the 
infinitesimal generator of the optimal wealth process in the Merton problem. Then, under the 
Inada conditions satisfied by the utility function U, we expect that this process does not reach 
the origin. Hence, we only need appropriate growth condition to ensure uniqueness. fDJ 

Remark 3.3 The first corrector equation has the following stochastic representation as the dy- 
namic programming equation of an ergodic control problem. For this representation we fix (i, z) 
and let {M','\t > 0} be non-decreasing control processes, for each i,j = Q,...,d. Let p be the 
controlled process defined by, 

pl = pl+Y^a}^'(s,z)Bi + Y^{Mr^Ml^^), 

./=i j=o 

for some arbitrary initial condition po and a d dimensional standard Brownian motion B. Then, 
the ergodic control problem is 

d{s,z) := inf J{s,z,M), 



where 



M 

\ [I l-T J d 

J(s,z,M) :=limsup - E - / \a'^(s)p,\ dt+ T A''^M^' 
7-^^ r L2jo 



i,j=0 

In the scalar case, this problem is closely related to the classical finite fuel problem introduced 
by Benes, Shepp & Withenhaussen p^. We refer to the paper by Menaldi, Robin and Taksar 
ll27l for the present multidimensional setting. 

The function w is the so-called potential function in ergodic control. We refer the reader 
to the book and the manuscript of Borkar [Tl [81 for information on the dynamic programming 
approach for the ergodic control problems. fDJ 



Remark 3.4 The calculation leading to ( 13.71 ) is used several times in the paper. Therefore, for 
future reference, we summarize it once again. Let v, 6, z and £, be as above. For any smooth 
functions 

: K'^ X M+ K, CJ ; R'j X M+ X W' ^ K, 

and e e (0, 1] set 

l"'is,x,y) ■=v{s,z)-e^(l)is,z)-e^(D{s,z,^). 
In the above calculations, we obtained an expansion for the second order nonlinear operator 

= e2(^-^|(jT^|2 + iTr[aaTc7^^]-^(/)+^^), (3.13) 

where a, £/ are as before and ^^{s,x,y) is the remainder term. Moreover, is locally 
bounded by a £ times a constant depending only on the values of the Merton function v, (j) and 
UJ. Indeed, a more detailed description and an estimate will be proved in one space dimension 
in Section|6l ITU 
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4 Corrector Equation in one dimension 



In this section, we solve the first corrector equation explicitly in the one-dimensional case. 
Then, we provide some estimates for the remainder introduced in Remark l334l 

4.1 Closed-form solution of the first corrector equation 

Recall that w = T7v,m>, a = tjvm, and the solution of the corrector equations is a pair (vv,a) 
satisfying, 

max{-^a^p^-^a"wpp+a,-V''^ + wp,~X°^^~wp}=0, w{s,z,0)=0, (4.1) 

where a = a/rj and oc(s,z) is given in ( I3.8l i. We also recall that the variables {s,z) are fixed 
parameters in this equation. Therefore, throughout this section, we suppress the dependences 
of (7, a and w on these variables. 

In order to compute the solution explicitely in terms of 77, we postulate a solution of the 
form 

( k4p^ + k2P^-+kip; Pi<p<Po 
w{p) - <^ vv(pi)-AO^i(p-pi); P<Pi (4.2) 
[ vi>(po) + A^^°(p-po); p>Po- 

We first determine k4 and ^2 by imposing that the fourth order polynomial solves the second 
order equation in (po,pi). A direct calculation yields, 

^4 = T7FT and ^2 = r^- 

We now impose the smooth pasting condition, namely assume that w is at the points po and 
p\. Then, the continuity of the second derivatives yield, 

, , 2fl _ /2a\i/2 

Po = Pi = ^ implymg that a > and po = -pi = (^^J . (4.3) 

The continuity of the first derivatives of w at the points po and p\ yield, 

4fc4(po)^+2A:2Po + ^:i = -A'''\ 

4fc4(pi)^+2fc2pi+^:i = A'-O. 

Since po = — Pi, we determine the value of ^1 by summing the two equations, 

Aio-Ao.i 

k\ ~ . 

' 2 

Finally, we obtain the value of a by further substituting the values of ki„ k^ and po = — Pi. The 
result is 



-a^^^pl and p^ = (g(A''« + A«^')) (4.4) 

All coefficients of our candidate are now uniquely determined. Moreover, we verify that the 
gradient constraint 

-A^^''<H'p <A'''i (4.5) 

holds true for all p e M. Hence, vv constructed above is a solution of the corrector equation. One 
may also prove that it is the unique solution. However, in the subsequent analysis we simply 
use the function m> defined in ( I4.2l i with the constants determined above. Therefore, we do not 
study the question of uniqueness of the corrector equation. 
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Remark 4.1 In the homothetic case with constant coefficients r, /i, and a, one can expUcitly 
calculate all the functions, see Section [8] Here we only report that, in that case, all functions are 
independent of the 5— variable and po,fl(z) are constants. Therefore, a{z) is a positive constant 
times the Merton value function. nu 

Remark 4.2 Pointwise estimates on the derivatives of w will be used in the subsequent sections. 
So we record them here for future references. Indeed, by (I4.5l l and the fact that w{-,Q) = 0, 

\wis,z,^)\<Iv,{s,zM\, \w^is,z,^)\<Iv,{s,z), where Ii^A^'WA'-O. (4.6) 

Moreover, under the smoothness assumption on v, we obtain the following pointwise estimates 

(|w| + |w,| + |w,,| + |w,| + \w,,\){z,^) < C{s,z){l + 1^ I), (4.7) 
I + 1^-.^ + 1^^-^ ^ Cis,z) and |w^^ | < (Cl[^„.^,]) (4.8) 

where C is an appropriate continuous function in M^, depending on the Merton value function 
and its derivatives. FDJ 



4.2 Remainder Estimate 



In this subsection, we estimate the remainder term in Remark [3741 So, let ^"^ be as in Remark 
We have seen in ( 13.131 1 that 

where a, £/ are defined in ( I3.8b - (l3.9t . and is the remainder By a direct (tedious) calcula- 
tion, the remainder term can be obtained explicitly. In view of our previous bounds (I4.7b - (l4.8l l 
on the derivatives of w, we obtain the estimate, 

\M'{s,z,^)\ < -e[^{^~r)m+^-aHee + 2^e)\c^^\]{s,z) 
+eC{s,z){l + e\^\ + e^\^\^ + e'\^\'), 

for some continuous function C{s,z)- Since U is C' and convex, 

\^'{s,z,^)\ < -e[^i^-rM\+^-aHe^^ + 2^eU,\){s,z) 
+eC{s,z){l+e\^\ + e^\^\He'\^\'), 

+{\(p-\ + e^\(p,\ + e0,|CJ^ |)|f7'(v- + £2|0,| + e4|c7-| + £3 0-|£i7^ |) - u'{v,)\ 
Suppose that GJ satisfies the same estimates (I4.7b -( l4.8b as w. Then, 

+eC{s,z){l+e\^\ + e^\^\^ + e'\^\'), 

+e^{\(j),\ + eC{s,z){l+e\^\)fu"{v, + e^\(j),\ + e^Cis,z)il+e\^\)). 
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5 Assumptions 



The main objective of this paper is to characterize the limit of the following sequence, 

Our proof follows the general methodology developed by Barles & Perthame in the context 
of viscosity solutions. Hence, we first define relaxed semi-limits by, 

«*(C):= limsup «^(C'), «*(C):= liminf S'(C') 

(£,?')^(0.f) (£.?')^(0,f) 

Then, we show under appropriate conditions that they are viscosity sub-solution and super- 
solution, respectively, of the second corrector equation ( 13.121 ). 
We shall now formulate some conditions which guarantee that 

i. the relaxed semi-limits are finite, 

ii. the second corrector equation (13.12b verifies comparison for viscosity solutions. 

We may then conclude that u* < m*. Since the opposite inequality is obvious, this shows that 
M = M*=M* is the unique solution of the second corrector equation ( 13.12b . 

In this short subsection, for the convenience of the reader, we collect all the assumptions 
needed for the convergence proof, including the ones that were already used. 

We first focus on the finiteness of the relaxed semi-limits and u*. A local lower bound 
is easy to obtain in view of the obvious inequality v^{s,x,y) < v{s,x+y) which implies that 

> 0. Our first assumption complements this with a local upper bound. 

Assumption 5.1 (Uniform Local Bound) The family of functions is locally uniformly bounded 
from above. 

The above assumption states that for any Vo := {s(i,XQ,yQ) e M_|_ x with xq +yo > 0, there 
exist ro = ro(vo) > and Eq = eo(vo) > so that 

fe(Vo) :=sup{ M''(i,A:,y) : £ Bro(vo), £ £ (0,eo] } < (5.1) 

This assumption is verified in Section |7] under some conditions on v and its derivatives by 
constructing an appropriate sub-solution to the dynamic programming equation ( 12.2b . However, 
the sub-solution does not need to have the exact behavior as needed in other approaches to 
this problem starting from 1361 l20l . Indeed, in these earlier approaches, both the sub and the 
super-solution must be sharp enough to have the exact limiting behavior in the leading term. 
For the above estimate, however, this term needs to be only locally bounded. 

The next assumption is a regularity condition on the Merton problem. 

Assumption 5.2 (Smoothness) The Merton value function v and the Merton optimal invest- 
ment strategy 9 are twice continuously dijferentiable in the open domain (0,°°)^ andv^{s,z) > 
for all s,z>Q. Moreover, there exists cq > such that 

CO < < 1 -co, Vi,zeR+. (5.2) 
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In particular, the above assumption implies that the diffusion coefficient a{s,z) in the first 
corrector equation is non-degenerate away from the origin. For later use we record that there 
exist two constants < a* < a* so that 

0< a* < -^^^ < a*, Vi,zGR+. (5.3) 
z 

We will not attempt to verify the above hypothesis. However, in the power utility case, the 
value function is always smooth and the condition (15.21 ) can be directly checked as the optimal 
investment policy 6 is explicitly available. 

We next formulate a natural assumption which was verified in ll36l in the context of CRRA 
type utiUties. This assumption will be used for the proof of the sub-solution property. To state 
this assumption, we first introduce the no-transaction region defined by, 

J^' ■.= {{s,x,y)eKe ■.Ali-Dv'{s,x,y)>0,andAlQ-Dv'{s,x,y)>Q}. (5.4) 

By the dynamic programming equation ( 12.2b . the value function v*^ is a viscosity solution of 

J3v^-^v^-t7(v^)=0 on 

Assumption 5.3 (No transaction region) The no-transaction region ,yV^ contains the Merton 
line ■.^{{s,z-e{z),e{z)) : s,zeR+}. 

Finally we assume that the second corrector equation (13.12b has comparison. Recall the 
function u introduced in (11.21 ). let b be as in ( 15.1b . and set 

B{s,z):=b{s,z~e{z),eiz)), s,zeU+. (5.5) 

Assumption 5.4 (Comparison) For any upper-semicontinuous (resp. lower-semicontinuous) 
viscosity sub-solution {resp. super-solution) u\ {resp. U2) of (13.12b in {0,°°)-^ satisfying the 
growth condition |m,| < B on (0,°°)^, i = 1,2, we have ui < u < U2 in (0,°°)^. 

In the above comparison, we do not make assumptions on the value of these functions at 
the origin. As discussed earlier, the operator £/ is the infinitesimal generator of the optimal 
wealth process in the limiting Merton problem. So we implicitly assume that this process does 
not reach the origin with probability one. 

6 Convergence in one dimension 

For the convergence proof, we introduce the following "corrected" version of H^, 

u^{s,x,y) := u^{s,x,y) - e^w{s,z,^), s>0, {x,y)<=Ki;. 

Notice that both families and u^ have the same relaxed semi-limits u* and m*. 

Theorem 6.1 Under Assumptions \5.1\\5^\5.3\ and \5.4\ the sequence {m^}£>o converges locally 
uniformly to the function u defined in ( 11.2b . 

Proof. In the next subsections, we will show that, the semi-limits and u* are viscosity 
super-solution and sub-solution, respectively, of ( 13. 12b . Then, by the comparison Assumption 
15.41 we conclude that m* < m < m*- Since the opposite inequality is obvious, this implies that 
u* = M* = u. The local uniform convergence follows immediately from this and the definitions. 
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6.1 First properties 



In this subsection, we only use the assumptions on the smoothness of the limiting Merton prob- 
lem and the local boundedness of {u^}£. We first recall that 

A := X'-'yV-'. 

Lemma 6.1 (i) For all e,s > 0, {x,y) E Kg, u^{s,x,y) > —eXv^{s,z)\y — 9(s,z)\. In particular, 
u*>0. 

(ii) If in addition Assumption 15. IM iolds true, then 

Q <Ui,{s,x,y) <u''{s^x,y) <°° for all i,^,^ > 0. 

Proof. Since (ii) is a direct consequence, we focus on (i). From the obvious inequality v*^ {s^x^y) < 
v{s,x + y), it follows that u^{s,x,y) > —e^w{s,z,^), so that the required result follows from the 
bound (I4.5l l on together with w(-,0) =0. lUJ 

We next show that the relaxed semi-limits u* and depend on the pair {x,y) only through 
the aggregate variable z = x + y. 

Lemma 6.2 Let Assumptions \5.1\ and \5^ hold true. Then, u* and are functions of (s^z) only. 
Moreover, for all s,z>0, 

u,is,z)= liminf M^(i',z' - 

and 

u*{s,z)= limsup m'^(/,z'-0(z'),0(z'))- 

(£j'.z')^iO,s.z) 

Proof. This result is a consequence of the gradient constraints in the dynamic programming 
equation ( 12.21 ). 

Af • (v^,v^) > and Ag j ■ (vf ,Vy) > in the viscosity sense. 
1. We change variables and use the above inequalities to obtain 

(l+A' V(l-a))v| >-A''Vvf, (l+A'^'£3a)v^ <A'''>eS'!, (6.1) 

in the viscosity sense. Since v*^ is concave in {x,y), the partial gradients and exist almost 
everywhere. By the smoothness of the Merton optimal investment strategy 6, this imphes that 
the partial gradient vf also exists almost everywhere. Then, by the definition of u^, we conclude 
that the partial gradients ii? and m| exist almost everywhere. In view of Condition (15.2b in 
As sumption 152] we conclude from ( 16.11 ) and the fact that vf > that 

<AeS'f. (6.2) 
We now claim that 

< f{s,x,y) 

:= v-{s,z - e) + e[u'^{s,x-e,y) + u^{s,x,y - e)) (6.3) 
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We postpone the poof of this claim to the next step and continue with the proof. Then, it follows 
from (16.2b . ( 16.31 ) together with Assumption l5.2l and ( 14.51 ). 



\ul{s,z,^)\ < e2A(v,(,,z) + v!(^,z,^)) 

< e^~Xiv,{s,z) + fis,zA))- (6.4) 

Hence, 

(ei -eo) = < el {v,{s,z) + f{s,z,^)). 

By the local boundedness of {w'^je, for any {s,x,y), there is an open neighborhood o/K of {s,x,y) 
and a constant K, both independent of £, such that the maps 

t u'^{s,x-t,y + t) + eKt and t ^ -u'^{s,x-t,y + t) + eKt 

are nondecreasing for all e > 0. Then, it follows from the definition of the relaxed semi-limits 
that u* and m* are independent of the (^-variable. 

2. We now prove ( 16.3b . For e > and {x,y) , {x ~ e) , {x,y — e) G K^, we denote as usual z=x+y 
and ^ = {y — 9{s,z))/£- By the concavity of in the pair {x,y) and the concavity of the Merton 
function v in z that: 

Vv(s,.^,3') < ^{v^{s,x,y)-v^{s,x-e,y)) 

< ^ {v{s,z) - v{s,z - e)) + ^ - e) - v^{s,x - e,y)) 

< v,is,z - e) + ^ {K^^z - e) - v%s,x - e,y)). 

By the definition of u^, 

v^^{s,x,y) < v,{s,z - e) + e{u^{s,x - e,y) + e^w{s,z- E,^e)) 

where (y- e))/e = ^ + - e))/e. We use the bound (14.6b on w, 

to arrive at, 

v?(5,;.,3') < v,{s,z -e) + eu'is,x -e,y) + e'XvMz){l + 1^1 + llM^^^llZlA^ . 
By exactly the same argument, we also conclude that 

v;is,x,y) < v,is,z^ e) + eu%s,x,y ~ e) + e'Iv,{s,z) (l + 1^ | + I ' ^ + ^(^'^)^- j . 
Then, using the bounds on 0- from Assumption |5.2| 

y!(.,z,^) - d,v'{s,z-e^-eis,z),e^ + e{s,z)) 

= (1 - e,is,z))v'^is,x,y) + e-Ss,zX.is,x,y) 
< Vz{s,z- e) + e{u'^ {s,x - e,y) + u'^{s,x,y - e)) 

+e^I.,is,z) (l + |a(.,z)| + 1^ I + l^i^^^)-^}^-'^) 



3. The final statement in the lemma follows from ( 16.4b . the expression of 7^ in (16.3b . and 
As sumption 15. II mj 
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6.2 Viscosity sub-solution property 

In this section, we prove 

Proposition 6.1 Under Assumptions \5.1\ and \5.2\ the function u* is a viscosity sub-solution of 
the second corrector equation ( 13.12b . 

Proof Let (io,Zo, <P) e (0,°°)^ x C^(M^) be such that 

[u* -<p){sq,zq)> {u* -<p){s,z) for all s,z>Q, {s,z)i^{sq,zo). (6.5) 

Our objective in the following steps is to prove that 

s^<p{sQ,ZQ)-a{so,ZQ) < 0. (6.6) 

1. By the definition of u* and Lemma l6^ there exists a sequence {s^ ,z'^) so that 

(/jZ*^) — > {so,zo) and u^{s'^,z^,0) — >u*{so,zq), as e\0, 
where we used the notation i3.2i . Then, it is clear that 

^^=^2^(/,^^0)-(p(/,z^) (6.7) 

and 

{x',/) = {z'~e{s\z'),0{s',z')) — > {xo,yo) {zo- e{so,zo),e{so,zo)). 

Since (u^) is locally bounded from above (Assumption 15.1b . there are ro := ro{sQ,xo,yQ) > 
and go := eo(so,^o,3'o) > so that 

b^:^sup{u^{s,x,y) : {s,x,y) e Bo, e e {0,eo]} <°°, where Bq := Brg{so,xo,yo) (6.8) 

is the open ball centered at {so,xo,yo) with radius tq. We may choose ro < zo/2 so that Bq does 
not intersect the line z~Q- For e,5 E (0, 1], set 

V'i^.z,^) := v{s,z)~e^£l^e^(p{s,z)-e\l + 5)w{s,z,^)~e^^'{s,z,^), 

where, following our standard notation (13.2b . 0*^ is determined from the function, 

fis,x,y) := C[is~s')' + {x + y~z')' + {y-eis,x + y))''], 

and C > is a large constant that is chosen so that for all sufficiently small e > 0, 

<j)^ > I +bt-(p, on Bq\Bi with Bi Br„/2(io,Jco,yo)- (6.9) 

The constant C chosen above may depend on many things including the test function (p,so,zo,S, 
but not on e. The convergence of {s^,z^) to {so,zo) determines how small e should be for ( 16. 9b 
to hold. 

2. We first show that, for all sufficiently small e > 0, 5 > 0, the difference {v'^ — or 
equivalently, 

re.S( ^ v^is,x,y)~\j/'''^{s,x,y) 
1 [s,x,y) .-— ^2 

= -u''{s,x,y) + (p{s,z)+£l + (j)''is,x,y) + e^5w{s,z,^), 
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has a local minimizer in Bq. Indeed, by the definition of u^, and i^, ( 16. 9t . (16.8b . and the 
fact that w > that, for any {s,x,y) £ dBo, 

I''^{s,x,y) > -u'{s,x,y)+il + l+b, + e^5w{s,z,^) >l+el> 0, 

for sufficiently small £ in view of ( 16.7b . Since I'^'^ (s^ ,y^) = 0, we conclude that 7*^'^ has a 
local minimizer {S^,x^ ,y^) in Bq with := x*^ +y^, := (y*^ — 0(s^,z^))/e satisfying, 

min (v^-r'^) = (v^-r'^)fe,|,)<0, |,~^-,o| + |z'^-zo|<ro, |^.|<n/e, 

(s,z.^)GBl 

for some constant ri . Since v*^ is a viscosity super-solution of the dynamic programming equa- 
tion ( 12.2b . we conclude that 

(^pv'-^W''^-U{w''^)yf,x',f) > 0, (6.10) 

and 

Ato.(v/^•^v^■^)(^^x^f) ^ (vA^■^-(l-A>V)v/f■^)(^^i^f) >o, 

By direct calculation using the boundedness of {s^ ,e£,^), we rewrite the last gradient in- 
equalities into the following 

-4e2(e|'^)3 + e^'z(/,f)[A'^°-(l + 5)vvp(/,^^p^)]+o(e3) > 0, (6.11) 
4e2(e|^)3 + eV(/,r)[A°^' + (l + 5)Wp(i^,f ,p'^)]+o(e3) > o, (6.12) 

where p*^ 1*^/7] (i*,?*^). 

3. Let po(i,z) be as in (14.3b . In this step, we show that 

\p^\<PQ{f,f) for all sufficiently small EG (0,1]. (6.13) 

Indeed, assume that p*^" < —po{s^",z^") = pi {s^" ,z^") for some sequence e„ G (0, 1] with e„ ^ 0. 
Then, Wp{s^" ,z^" ,p'^") = — A^ ^ and it follows from inequality ( 16.121 ). together with the fact 
P*^" < piiS^",z^") < 0, that 

< 4e2(e„|'^")3 - eSis'",f")5X^^-^ + o{e^,) < -En^v,if",f")5X^^-^ + o(e„3). 

Since 5 > Q, this can not happen for large n. Similarly, if p'^" > Pq{s^",z^" ) for some sequence 
£„ 0, we have Wp{S^" ,z^" ,p'^") = A' and it follows from inequality ( 16.1 lb . together with the 
fact that p^" > po(i'^",?") > 0, that 

0< -4£2(£„|^")3 + e>,(/",f^")(-5A'''') + o(e„3) < -eSiS"" ,f")5V'' + o{e'„), 

which leads again to a contradiction for large n, completing the proof of ( 16.131 ). 

4. Since {s^,z'^J is bounded and {s,z) H> Po{s,z) is continuous, we conclude from (16.13b that 
the sequence {^^)e is bounded. Hence, there exists a sequence e„ — > so that 

(i„,z,„^„) f'Sl^'O — > (i,£,|)==(io,zo,|) 

for some ^ G M. The fact that the limit of {s„,z„) is equal to {so,zq) follows from standard 
arguments using the strict minimum property of {sq ,zo) in (16.5b . We now take the limit in (16.10b 
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along the sequence (e„). Since the function has the form as in Remark [3741 we do not 
repeat the computations given in Section[3]and, given the remainder estimate of section l4!2l we 
directly conclude that 

£„->0 \ / 

= ^{ri<7^){sQ,ZQ)i^ + ^{l + 8)a-(so,zo)w^^{so,zo,^)~£/(p{so,zo) (6.14) 

In the above, we also used the fact that all derivatives of vanish at the origin as e tends to 
zero. 

5. In the Step 3, we have proved that \p£\< Poize). Hence, |^| < (i?Po)('50;-Zo)- Since w — r]VzW, 
a = rjv-a, the first corrector equation ( 13.1 11 1 implies that 



a{so,ZQ) = ^i<y^ri){so,zo)^^ + ^a^{so,zo)w^^{so,zo,^)■ 



'We use the above identity in ( 16.141 ). The result is 



,zo) < ^{a^Tl){so,zo)l,^ + ^{l + 5)a^{so,zo)w^^{so,zoA) 
1 



= a{sQ,zo) + ^5a^{so,zo)w^^{sQ,zo,^)- 

Finally, we let 5 go to zero. However, | = | ^ depends on 5 and care must be taken. But since 
l^nl < {'nPo){smZn), it follows that |^ is uniformly bounded in 5. Hence the second term in the 
above equation goes to zero with 5, and we obtain the desired inequality ( 16.6b . FDJ 

6.3 Viscosity super-solution property 

In this section, we prove 



Proposition 6.2 Let Assumptions pJ] 15.21 and \5.3\ hold true. Then, the function m* is a viscosity 
super-solution of the second corrector equation ( 13.12b . 

We first need the following consequence of Assumption ^ .3 1 and the convexity of v"^. Similar 
arguments are also used in 



Lemma 6.3 Assume the hypothesis of Proposition |6!2l Let (jc,y) be an arbitrary element of K^. 
Then, 

(i) fory > 9{s,z) (or equivalently, B, > 0), we have Aq j ■ {v^{s ,x,y) ,Vy {s ,x,y)) > 0, 
(u)fory < 9{s,z) (or equivalently, B, < 0), we have Aj q ■ {v^{s,x,y),Vy{s,x,y)) > 0. 

Proof. For z G IR.+ set 

elis, z):= sup {y:{z-y,y)eK,, and Ag i • (v^, vp(i,z-y,3;) = O}. 

In view of the form of Kg, we have y > — z/(e^A'''' ) and by convention the above supremum is 
equal to this lower bound if the set is empty. By the concavity of v*^, we conclude that 

A^ -(v' v'](, xv]^ =Oforally< 0^(i,z), 
Ao,i [v^,Vy)[s,x,y)< >oforally>0^(5,z). 
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Let ,yy^ be as in (15.4b . Therefore it is included in the set {{s,x,y) : y > 0^{s,z)}- Since As- 
sumption lS . 3 1 states that the Merton line {{s,x,y) : y = 9{s,z)} is included in y^^, we conclude 
that 9{s,z) > 9^{s,z)- This proves the statement (i). The other assertion is proved similarly, nu 

Proof of Proposition^^ Let {sd.zq,^) G x be such that 

{ii^-(P){sq,zq) <{u^.-(p){s,z) for all s,z>0, (s,z) ^ (io,Zo)- (6.15) 
Our objective in the following steps is to prove that 

£^(p{sQ,Zi))~a{sQ,ZQ) > 0. (6.16) 

1. By the definition of u^, and Lemma |6^ there exists a sequence {s^ ,z^) so that 

(/jZ*^) — > {so,zo) and m'^(/,z'^,0) — >u^,{so,zq), as e\0, 
where we used the notation i3.2i . Then, it is clear that 

4 :=^2^(/,^^0) -(?(/, z^)^0 

and 

^ {z' - 0(/,z'^)) ^ {xo,yo) (zo - 0(io,zo), 0(^o,zo)). 

Since u^{s,x,y) > —e^w{s,z,^) > —£C{s,z)\y — 9{s,z)\, for some continuous function C, there 
are tq ro{so,xo,yo) > and £0(^0,^0, Jo) > so that 

b* -.^ inf M'^(5,x,y) > -00, where Bq := Brg{so,xo,yo)- 

{s.x,y)eBo 

We also choose ro sufficiently small so that Bq does not intersect the line z = 0. For e e (0, 1] 
and 5 > 0, define 

r-\^,z,^) := v{s,z)-£^i: - e2(p(5,z) -£4(1 - 5)w{s,z,^) + £^'(s,z,^), 

where, following our notation convention ( 13.2b . the function is obtained from the function 
(j)^ defined by, 

^'{s,x,y):=C[{s-s')U{x + y-z')' + {y-9{s,x + y)f]] 
and, similar to the proof of the super-solution property, C > is a constant chosen so that, 

-b* + ei+{(p-(j)''){s,x,y) <Q on dBo- (6.17) 

2. Set 

I'-H^.z,^) ■■= e-\v'-xt,'-')is,x,y) 

= -u'^{s,x,y) + (p{s,z) + il - ^^is,x,y) - e^dw{s,z,^). 

Since w{s,z,0) = 0, we have I^'^ {s^ ,z^ ,0) ~ 0. On the other hand, it follows from ( 16.17b that 

l'=-^{s,z,^)<-b* + £l+{(p-(p')is,x,y)<-e'5w{s,z,£,)<0 on dBo. 
Then, the difference v*^ — 1//*^^^ has an interior maximizer (i*,!*^, 1*^) in Bq, 

max(v^ - xir^-') = (v^ - y/^-'){f ,xe,f), and -io| + |z£ -zo| + He\ < n, (6.18) 
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(6.19) 



for some constant ri. By the sub-solution property of v"^, at {S^,X£,y£), 

3. In this step, we show that for all sufficiently small e > 0, 

Ag,l•(V/f^V<•^)(^",i^,/.)>0 and Al^- {w','^Wy-^){f ,^~eJe) > 0. (6.20) 
By Lemma |631 it suffices to prove that 

DO•l:=Ag l•(v/,''^r^■^)(^^i■„/,)>0 for | < 0, 
D'O Af „ . (v/f^ > for | > 0. 

We directly compute that 



(6.21) 



V;-e>;-e4(l-5)(w,-^vv^ ' '''■2^^^' -^^^ ^ ^^^^ 



e 

1-a 



^) +4e^C{{z-z')^ + (1 - 0,)(.v- 0)^ 



Then, it follows from the estimates (16.18b that 



,0,1 



e^( + (1 - 5)vv^ +A''''v,)(/,lM^) -4Ce2(e|^)3 + o(e3) 
Z)''° = e^(-(l-5)vv^+A''V)(/,zM^)+4Ce2(e|^)3 + o(e3). 

Since w solves (|4l]i, +A''''v; > and -w^ + A' °V; > 0. Then, 



,0,1 



and 



D 



> -e^dw^if,f,^'') + o{e^) for | < 0, 

1'" > e35w^(/,zM^)+4Ce2(£|^)3 + o(£3). 

> e^5w^(5^,zM^) + o(e3) for | > 0. 



Hence, ( 16.21b holds for all sufficiently small e > 0. 

4. In this step, we prove that £,e is bounded in e e (0, 1]. Indeed, in view of ( 16. 19b and ( 16.20b . 
we conclude that 



> [l5v'-^Y''^^U{xi/!'^)){s',x„y,) 



\^e\^ + K^aHs',Zs)w^^{zs,ls) 



2 ' 2 

{s^,Ze)+^'{S',x„y,) 



(6.22) 



where we used the fact that the function xj/^^^ is exactly as in the form assumed in Remark [3!4l 
Then, by the remainder estimate of section l4!2l we deduce that, 



''{s',x„ys)\<C{s',Ze) e + e|^.| + e^|^, 



(6.23) 



In Section m the function w is explicitly constructed. Since w is linear in ^ for large values of 
^, there is a continuous function C{s,z) so that 



0<w^^{s,z,^)<C{s,z), for all £ 
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Then, since (S^jZe) is uniformly bounded in e G (0, 1], there are constants C,C > so that: 



0>£^C 



2|S |2 



Hence {£,£)£ is also uniformly bounded in e G (0, 1] by a constant depending only on 5 and the 
test functions. 

5. Since {ze,^£)E£[0A] is bounded, there exists a sequence (e„)„ such that 

£„ \0 and (z,„^„) (z,„,^,„) (z,|) = (zo,|) £ (0,-) x M, 

where z = zo follows from the strict maximum property in ( 16.15b by a classical argument in the 
theory of viscosity solutions. We finally conclude from (16.22b and (16.23b that 

> ~^ia\^){so,ZQ)^^-£/(pisQ,ZQ)-£/^iO) + ^(1 - d)a^{so,ZQ)w^^{so,zo,^) 
= -£/(p{so,zo)-'^{<y^v,,,){so,ZQ)i^ + ^{l-d)a^{so,zo)w^^{so,ZQ,^), 
since (0) = 0. Now, from the first corrector equation ( 13.111 ), this implies that: 

> ~£/(p{so,zo)+a{so,zo) + ^8a^{sQ,zo)w^^{so,zo,£,)- 
Finally, we conclude that £/(p{so,zo) — a{so,zo) > 0, by sending 5 to zero. fDJ 



7 Verifying Assumption I5J 



In this section, we verify Assumption (15.1b . This is done by constructing an appropriate sub- 
solution of the dynamic programming equation ( 12.2b . Clearly, this construction requires as- 
sumptions and here we present only one possible set of assumptions. To simplify the presenta- 
tion, we first suppose that the coefficients are independent of the i-variable. Next, we assume 
that there exist constants < A;* < A:* so that the limit Merton value function satisfies 

0<A*z< T7(z) <rz. (7.1) 

Let c be the optimal Merton consumption policy given as in (12.6b . We assume that 

U{c{z))>k,zv'{z), (7.2) 

for some constant > 0. Notice that all the above assumptions hold in the power utility case. 
First, using ( 15.3b and the explicit representation of a, one may directly verify that there is a 
constant a* > so that 

a{z) < a*zv'{z). 

Then, the definition of £/ and the above assumptions imply that 

^y(z) = U{c{z)) > k,zv'{z) > ^fl(z) = '^^u{z). (7.3) 
In view of the comparison assumption, we conclude that 

0<m(z) < ^ v(z). (7.4) 
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Moreover, since we assume that coefficients are independent of the s variable, (12.71) is equivalent 
to 9{z) = 77(z)(jLi -r)/(j2. Hence, i5.3i implies that 

-v"(z)<T7(z)v"'<-2v"(z). (7.5) 

We now use these observations to construct a sub-solution of the dynamic programming 
equation of the form 

V'{x,y) := v(z) - Ke^iz) + £'W{zA), (7.6) 
with a sufficiently large constant K > a* /k^ and a slightly modified corrector, 

W(z,^):=zv'(z)w(^/z), 

where the function w{z) and the constant a > are the unique solution of w(0) = and 

_^p2_ (^^^^+5- „2A'-« + Wp; -2A«'i-Wpl. (7.7) 



The solution of the above equation is explicitly available through the general solution obtained 
earlier in Section [43] 

The fact that V'^ is a sub-solution of ( 12.21 ) follows from a tedious but otherwise direct calcu- 
lation. To streamline these calculations, we first state an estimate that follows from the explicit 
form of W. 

Lemma 7.1 There is a constant k* > so that 

zm^iz,^)\<kW{z), 
|W,(z,^)|<Fv'(z)(l + f), 

z\d,W{z,^)\ +z\dyW{z,^)\ < k*zv'{z) (i + f ) ' 
Proof. These estimates follow directly from straightforward differentiation and the estimates 

dLB-dlSl)- nu 



Lemma 7.2 (Lower Bound) Assume (17.1b , ( 15.21 ) and i2.6\ . Then, for sufficiently large K > 0, 
defined in i7.6i is a sub-solution of ( 12.21 ) in M^. Moreover, 

u'ix,y)<Kv{z)+£^Wiz,^) 
on and Assumption [5J] holds. 



Proof. We need to show that at any point {x,y) G one of the three terms in (I2.2l i is non- 
positive. Since {x,y) E R^, by assumption (15.2b . we have 

|^H^^<^. ^ s:^iel[-i,i]. 
e e z £ 



Let po > be the threshold in the equation (17.7b . We analyze several cases separately. 
Case 1. Po < S < 1/g. 
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In this case, (z, ^) = 2X^'^v'{z). We use the previous Lemma and (15.2b . to arrive at, 



4.0 



1 



e 



.3-1 1.0.,/ / 1 , ,*„3\ 



< e'V%'{-l+k*e')<0, 
provided that e is sufficiently small. 

Case 2. -l/e<S<-po. 

A similar calculation, shows that Aq j • (y^,V^) < 0, for all sufficiently small e. 

Case 3. \Z\ < po. We now use Remark|33]to conclude that 



^ ^ + -f^ W^^ (z, ^ ) - K^v{z) +M'{z,^) 



2^2 

We first use (IXTT i. (EU, dTTI i, dTJl l and set p := ^ /z. The result is 



< e'v'(z)77(z) 



(a*F)2 



-e-^'^fz, 



= e^''(z)T7(z) +e2^^(z,^). 

If AT is sufficiently large then K{k^,)^ is larger than a and by (17.11 ). the above estimate impUes 
that 

y <-zv'iz)+^'iz,^). 

We now estimate by recalling the results of subsection 14.21 We split this in three terms 
coming from the value function v, the corrector W and from the utility function. 



We estimate each one using Lemma im Then, 



Mf. < K 



eZ{n r)zv'{z) + ^ H 2 + 2eS(0/z)) zy'{z) 



< eKk*zv'{z). 



Also 



PW -rz{{l-{e/z)) + eZ)W,-fiz{eZ + {e/z))Wy 



Finally 



< k*zv\z). 



n = U{v')-U{V,') 

< U{v')-U{v'[l~e^K + k*e^])<0. 
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Hence, there is k* so that. 

I^^l <erzv'(z). 

Hence if K is sufficiently large, is a sub-solution for all small e. 
Boundary y ~0. 

Then, again by (I5.2l i. for all sufficiently small e > 0, 

„ y-0{z) -d{z) 

^ = ^— = -^<-po. 

Hence, by the second case, and Lemma 163] 

Af,o • (v/,y/)(^,o) < = Af,o • (v^ vp(^,0), > 0. 

Boundary X = 0. 

By a similar analysis, we can show that 

A|;^i • iV,^,Vf)iO,y) < = Ai, ■ (v^ vP(0,y), Vy > 0. 

Then, on M^, is a sub-solution of i2.2i while is a solution. Also on the boundary of 
again is a sub-solution of an oblique Neumann condition and v*^ is a super-solution. Then, 
by comparison, we conclude that > on R^. This proves the lower bound on u^ on the 
positive orthtant. mj 

8 nomothetic case 

In this short section, we consider the classical CRRA utility function 

U{c):^-^, c>Q, 
for some y > with 7=1 corresponding to the logarithmic utility. Then, 



with the Merton constant 



Y 2 7' C7^ 



Hence the risk tolerance function and the optimal strategies are given by 

T7(z) = ^, 0(z) = — 2-z:=7rMZ, c(z) = vmZ- 
7 7C7^ 

Since the diffusion coefficient a{z) = CT0(z)(l — 0^-), 

The constants in the solution of the corrector equation are given by, 

1/3 



Po 
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Since 

£/v{z) = U{C{Z)) = (VMZ)'"^ = VMV(Z), 

the unique solution u{z) of the second corrector equation 

27 



is given by 



(t2(1-7) 

where 



u — /; 2 -1 / \ 1-r 
«(z) = 2y ^ "''^ "^^^ " ' 



Finally, we summarize the expansion result in the following. 
Lemma 8.1 For a power utility function U, 

v'{x,y) = v{z)-eh,oz'-^ + 0{e'). 
The width of the transaction region for the first correction equation — 2ri(z)pQ is given by 

The above formulae with X'--^ ~ 1 are exactly the same as equation (3.13) in ||20) . 
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